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Introduction 


In 2020, nations hastened to contain an emerging COVID-19 pandemic by deploying diverse public health approaches, 
but conclusive appraisals of the efficacy of these approaches are elusive in most cases. One of the medicines 
deployed, ivermectin (IVM), a macrocyclic lactone having biochemical activity against SARS-CoV-2 through 
competitive binding to its spike protein, has yielded mixed results in randomized clinical trials (RCTS) for COVID-19 
treatments. In Peru, an opportunity to track the efficacy of IVM with a close consideration of confounding factors 
was provided through data for excess deaths as correlated with IVM use in 2020, under semi-autonomous policies in 
its 25 states. 


Methods 


To evaluate possible IVM treatment effects, excess deaths as determined from Peruvian national health data were 
analyzed by state for ages >60 in Peru's 25 states. These data were compared with monthly summary data for excess 
deaths in Peru for the period 2020-2021 as published by the WHO in 2022. 'To identify potential confounding factors, 
Google mobility data, population densities, SARS-CoV-2 genetic variations, and seropositivity rates were also 
examined. 
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) IVM, the channel by which most patients had obtained it previously, were then stopped, with its further use allowed 
only by a doctor's prescription [8-12]. Nationwide changes in daily excess all-cause deaths before and after these 
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deaths both deposited and freely available in the Dryad data repository [13,14]. Notably, the February 23, 2021 data 
snapshot closely aligned with the most recently available Peruvian national health data and with WHO summary 
data for excess deaths in Peru, in the period 2020-2021, with the comparison procedures as detailed [14]. 


Conclusive assessments of the effects of drug distributions and treatments on a national scale, however, are typically 
precluded by confounding influences. One such potential pitfall with the analysis of population-level mortality data 
is the use of case fatality statistics, which are often unreliable and which were indeed underreported for COVID-19 in 
Peru [15]. In addition, the unstratified aggregation of all age groups could distort this analysis. Different states in Peru 
have varying age distributions, and the percentages of COVID-19 cases across age groups could change through the 
course of the pandemic. Another extraneous influence that could have affected COVID-19 mortality in Peru was a 
diversity of genetic lineages of SARS-CoV-2 that circulated at various incidence rates state by state in Peru and 
worldwide in 2020 [16,17]. Moreover, a social isolation mandate imposed nationally in May 2020 would have affected 
death rates from COVID-19 in Peru. Other confounding factors were varying seropositivity rates and population 
densities in the different Peruvian states, as detailed below in the Discussion section. 


The availability of robust data for these potentially confounding influences and for excess death values by age 
allowed an assessment of these distortions and the derivation of reliable conclusions. Moreover, facilitating analysis 
of correlations between IVM usage and excess deaths were robust variations in the extent and timing of IVM 
distributions in the 25 Peruvian states. In 10 of these states, IVM was distributed widely within brief periods that 
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states between April and August 2020. The typical dose of IVM provided to both COVID-19 inpatients and outpatients 
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treatments began there on a large scale [20]. 


Mega-Operaci6n Tayta (MOT) 

IVM was typically distributed through regional health offices, voluntary channels, and other private groups, as 
detailed for several states [20]. However, 10 states distributed IVM on a mass scale through a national program led by 
the Ministry of Defense, Mega-Operacidn Tayta (MOT). Two of these states had confounding factors for their 
distributions of IVM. Pasco had three different distribution dates, July 23, August 5, and August 25 [27-29], while 
Junin's MOT deployment, which began August 4, had been preceded by state distributions of IVM to health centers 
beginning July 22 [30,31], 13 days earlier. 


MOT, an extension of a precursor program, Operacidn Tayta [32], was spearheaded by the Peruvian Ministry of 
Defense and Army. Eleven other government agencies partnered in this effort, including the ministries of health, 
interior, agriculture, and education, while participating personnel included those from the Army, Navy, Air Force, and 
police [33]. MOT's objective was to reach every part of a targeted region using rapid response teams that partnered 
with local health officials. These teams detected COVID-19 cases house by house, administered IVM to patients and 
family members in their households, and gave them food to encourage their isolation for 15 days [34]. 
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2022 meta-analysis of IVM treatment studies for COVID-19 [62]. Remarkably, however, in contradiction to that 
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mediate the virus's morbidity [63-65]. This molecular mechanism could explain the sharp increases in 


SpO» (peripheral oxygen saturation) values observed in severe COVID-19 patients within 24 hours after IVM 


treatments in three studies, which tracked SpO, values without the use of supplemental oxygen [66-68]. 


IVM is practical for widespread distribution on a global scale, having been deployed in worldwide campaigns to 
eradicate two scourges, onchocerciasis and lymphatic filariasis [69]. It holds a 34-year record of safety [4,6,70], 
including at doses considerably higher than the standard dose of 200 ug/kg [71-73]. It was used in RCTs for COVID-19 
treatment at cumulative doses of 1,500 ug/kg [49], 3,000 ug/kg [74], and 6,000 ug/kg [75] over five days with only small 
percentages of mild or transient adverse effects. 


This article was published as a reprint on OSFPrints in 2021 [106]. The analysis presented here extends and refines 
the work reported in the preprint by the same authors. 


Materials & Methods 
Study design 


An ecological study design was used. The health tracking values used for the analysis were compiled daily by the 
Centro Nacional de Epidemiologia, Prevencidn y Control de Enfermedades (National Center for Epidemiology, 
Prevention and Disease Control) and Instituto Nacional de Salud (National Institute of Health, NIH) in Peru. The 


tracking value used for the analysis was deaths from all natural causes (excluding violent deaths), hereinafter 
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beginning in May 2020 fluctuated by more than double the baseline value in January through February 2020, 
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then also tracked at 30 and 45 days following the date of peak deaths. Figure 1 provides examples of graphs of seven- 
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official summarized, this lockdown “failed completely” because for 75% of Peruvian residents, “if they do not work 
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one day, they cannot eat” [2]. An indication of actual compliance to such official lockdown orders is provided by 
Google Community Mobility Report data from cell phones within a given locality, which allow an objective 
quantification of social interactions [78-81]. Actual versus mandated changes in social mobility elsewhere in the 
world were found to have likewise varied considerably during the 2020 pandemic period. In some countries, such as 
Sweden, certain mobility restrictions were undertaken on individual initiatives [79], while in others, official 
mandates had limited impact on actual mobility [80,81]. Therefore, to most objectively assess the potential effects of 
social isolation policies on mortality trends in Peru, six indices of the Google Community Mobility Report data were 
retrieved for each of Peru's 25 states and compared with mortality trends. 


Results 


Analysis on excess all-cause deaths was performed state by state, with the population ages >60. For each state, the 
date of peak excess deaths in its first wave of the pandemic, as specified in the Materials and Methods section, was 
determined. Decreases in excess deaths from the date of peak deaths to 30 and 45 days afterwards were then tracked. 
The 25 states of Peru were grouped by the extent of IVM distributions: maximal (mass IVM distributions through 
operation MOT, 10 states), medium (locally managed IVM distributions, 14 states), and minimal (Lima, with 
restrictive policies). Table 1 shows the mean reductions in excess deaths at 30 and 45 days after date of peak deaths 
by the tier of IVM distribution, as also shown graphically for +30 days after peak deaths in Figure 2C. Table 2 shows 
these values for each state individually. 
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Figure 2: Excess all-cause deaths, for the national population of Peru and by state, overview 

A) Excess all-cause deaths (all ages), national population of Peru. For A and B, the y values are seven-day 
moving averages; for B and C, ages >60. Data are from Peru's National Information System of Deaths 
(SINADEF) [76], with associated frozen datasets available from the Dryad data repository [13]. B) Reductions in 
excess deaths for all states of Mega-Operacidn Tayta (MOT), a program of mass ivermectin (IVM) 
distributions, excluding Pasco, which had them on three different dates. e MOT start date; A date of peak 
deaths; m 30 days after date of peak deaths. (Junin had a distribution of IVM 13 days before the MOT start.) C) 
Reductions in excess deaths at 30 days after peak deaths for the 25 states, by extent of IVM distributions: 
maximal- MOT (+), mean -74%; moderate-local distributions (O), mean -53%; and minimal-Lima (X), -25%. > 
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Note that for Lima, as detailed in a previous analysis [20], large distributions of IVM beginning in August marked the 
end of its prior policy of IVM restrictions, with a second peak in excess deaths occurring on August 4, followed by a 
decline, as shown in Figure 1. However, as per the methodology of this analysis, the date of Limas's first peak in excess 
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effect, as manifested by associated p-values much less than 0.05. These overall sharp correlations were obtained 
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was restricted support the possibility that IVM treatments for COVID-19 were effective in Peru in that timeframe. As 
shown in Figure ZA, between August 1 and December 1, 2020, during the period of IVM use, nationwide excess deaths 
decreased by 14-fold. After IVM use was restricted [13,75], however, excess deaths then increased 13-fold through 
February 1, 2021. 


Discussion 


The 25 states of Peru that conducted IVM treatments for COVID-19 at different time periods provide a robust set of 
subpopulations in which the treatment impacts can be evaluated. For the 10 MOT states, excess deaths dropped most 
sharply, by mean values of 74% at +30 days and 86% at +45 days after the date of peak deaths. For the 14 states with 
locally administered IVM distributions, excess deaths dropped by means of 53% at +30 days and 60% at +45 days. In 
Lima, however, where IVM treatments were delayed until August, four months after its initial pandemic surge in 
April, excess deaths dropped by only 25% at +30 days and by only 25% at +45 days after its date of peak deaths on May 
30. The MOT states had sharp drops in excess deaths after reaching peak values in close time conjunction after MOT 
start dates (Figure 2B). For the full set of 25 states, reductions in excess deaths correlated with the extent of IVM 
distribution, maximal, medium, or minimal, with p<0.002 using the (two-tailed) Kendall rb test. 


Given the association between IVM treatments and sharp mortality reductions revealed in this analysis, neither 
random fluctuation nor an unidentified, extraneous cause of these reductions in deaths appears likely. However, it is 
useful to consider the potential confounding influences of social isolation, changing seropositivity rates, variations 
in viral strains across states, and other factors. To begin with the most straightforward of these considerations, 
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could explain reduced a tă impacts there, several other IVM-treated states with low seropositivity rates had 
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[87]. As for people living in the same ia seal a demographic sava in 2017 showed that Lima households with 
more than five people represented 27% of the total; in Loreto, that figure was 42%, and in Ucayali, 36% [14,88]. Thus, 
neither population densities per area nor densities per household were markedly different in Lima versus population 


centers of other states for which this analysis was performed. 


It had been proposed that cross-immunity from the dengue virus, which causes dengue fever, could explain lower 
than expected levels of COVID-19 mortality in some regions of South America, including Brazil [89]. This theory 
collapses in Peru, however, with greater than 70% reductions at 30 days after peak deaths, for example, in both the 
Peruvian states of Moquegua, which has not had dengue cases in the last 20 years, and Loreto, the epicenter of 
dengue in Peru [90,91]. Finally, one other data artifact could be that several peaks and drops in Lima's different 
districts could explain the low reduction in excess deaths. However, the pattern of total deaths, ages >60, for most of 
these districts, those comprising the bulk of the population, is the same: rising deaths to a peak around late May 2020 
and then a three-month plateau following [14]. 


One potential confounding factor that cannot be conclusively resolved is the potential role in the operation MOT of 
the distribution of IVM for prophylaxis to members of the households of COVID-19 patients and others at high risk for 
the disease who had not tested positive for it (see section on MOTabove). A preventative effect of IVM against COVID- 
19 may have contributed to the sharp reductions in excess deaths that followed its use in the 10 MOT states. The 
potential efficacy of IVM for COVID-19 prevention was indicated in two human prophylaxis studies, both using IVM 
in doses of at least 150 ug/kg per week. These RCTs reported statistically significant (p<0.001) reductions of 87% and 
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On May 7, 2021, the date of this WHO report, the seven-day moving average of COVID-19 deaths per million of 
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population was roughly the same in Uttar Pradesh, all of India, and the United States, i.e., 1.4, 3.2, and 1.9, respectively 


Ivermectin Prophet Ugest fel IA 97 iDO7ilea Pisgesiipe7 Oberotatiasal Spui afftaș.128 Sieg insine Psfpenaitinseere Mateliloation of IVM, the seven- 


(articles/82162-ivergjă şir popi a avast 8 fe Epic dese, Atta ptiepa etuey 9 gipsi -U in SzpsepensițViscpiis-matchina) ati ve total of COVID-19 
Cp pl +7 more (/articles/82162-ivermectin-prophylaxiş-used-for-covid-19-a-citywide-prospective-observational-study-of-223128-subjects-using-propensity-score-matching? 
deaths per milion In population from July 7, 2U7I through April], was 4.3 în Uttar Pradesh, as compare with 

70.4 in all of India and 1,596.3 in the United States [100]. These data are from the Institute for Health Metrics and 


Lucy Kerr, Flavio A. 


authors-tab=true) 


idend9 9 IWPOVORESIA RATIONS/sIGN_ 


le) 


&agistrationsâdgutultă 


strictiy-controlled-păS stă $hedinited-Sţates in that period may reflect the use of these same home treatment kits containing IVM, 


Lucy Kerr, Fernandelop ci zcipreal/aatisits/ii1961i'moglonite othewmstites af Irvttiie-Ț99-404| Tisei piete A2-6f COX IDcoric tă nerietivertețeaf the seotber-menner-results-of-a- 


prospective obser 0 onaetus Dos But ret MAE ulațienrof-88012-subigqts rsrtab=tru in Uttar Pradesh and highest in the United States [100] 


cannot explain the trends in COVID-19 deaths that occurred. 


Unlike for Peru, where granular data by state and by time period for excess deaths allow a well-founded correlation 
of trends for excess deaths with the extent of IVM use, the data for Uttar Pradesh are suggestive but not conclusive 
for efficacy of IVM in the treatment and possibly also the prevention of COVID-19. The possibility of such efficacy 
having been achieved in Uttar Pradesh, however, is supported by the particularly high efficacy of IVM in conjunction 
with doxycycline and zinc in COVID-19 treatment (the combination used in Uttar Pradesh) as reported in two clinical 
studies [66,67]. 


The exceptional safety profile and low cost of IVM certainly support its use as in Peru's operation MOTand in Uttar 
Pradesh as an attractive national policy for COVID-19 mitigation. These significant reductions in mortality as 
achieved in Peru and Uttar Pradesh suggest that the impact of such a national IVM deployment would be observable 
within a relatively short period. However, generic drugs have often fared poorly in competition with patented 
offerings in past decades, based upon the unfortunate vulnerability of science to commodification and regulatory 
capture [102,103]. For example, the inexpensive triple therapy for peptic ulcers, proven 96% curative in 1990 and now 
the standard of care, was delayed for a decade in its worldwide deployment until patents for two best-selling 
palliative drugs for peptic ulcers expired [48]. Such a potential bias against IVM was suggested by a February 4, 2021 
press release from Merck, which was then developing its own patented COVID-19 therapeutic, claiming that there 
was “a concerning lack of safety data” for IVM [104]. However, IVM is Merck's own drug, found safe at doses 
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